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Abstract. In this paper we are interested in moments of Minkowski question mark func- 
tion It appears that, to certain extent, the results are analogous to the results obtained 
for objects associated with Maass wave forms: period functions, L-series, distributions, spec- 
tral properties. These objects can be naturally defined for T{x) as well. Various previous 
investigations of 7{x) are mainly motivated from the perspective of metric number theory, 
Hausdorff dimension, singularity and generalizations. In this work it is shown that analytic 
and spectral properties of various integral transforms of do reveal significant informa- 
tion about the question mark function. We prove asymptotic and structural results about 
the moments, calculate certain integrals involving ?(a;), define an associated zeta function, 
generating functions, Fourier series, and establish intrinsic relations among these objects. 

Mathematical subject classification: 11A55, 11M41, 11F99, 26A30. 
Keywords: Minkowski question mark function, period functions 

1. Introduction 

The aim of this paper to continue investigations on the moments of Minkowski F{x) 
function, begun in [Ij and [2j. The function F{x) ("the question mark function") was 
introduced by Minkowski in 1904 as an example of a continuous function F : [0, oo) — > [0, 1), 
which maps rationals to dyadic rationals and quadratic irrationals to non-dyadic rationals. 
For non-negative real x it is defined by the expression 

F([ao, ai, a2, ag, ...]) = 1 - 2-"» + 2-("«+"i) - 2-(«o+ai+a2) ^ 

where x = [00,01,02,03,...] stands for the representation of x by a (regular) continued 
fraction [13]. More often this function is investigated in the interval [0,1]; in this case we 
use a standard notation ?(a;) = 2F{x) for x G [0,1]. For rational x, the series terminates 
at the last nonzero element o„ of the continued fraction. This function was investigated by 
many authors. In particular, Denjoy [9] showed that ?(x) is singular, and the derivative 
vanishes almost everywhere. The proof of this fact is also given in The problem arose in 
connection with the Calkin- Wilf tree (see below); the Farey tree is actually a subtree of it. In 
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fact, singularity of follows from Khinchin's average value theorem on continued fractions 
([13], chapter III). The nature of singularity of F{x) was clarified by Viader, Paradfs and 
Bibiloni [22]. In particular, the existence of the derivative 7'{x) in R for fixed x forces it to 
vanish. Some other properties of ?(x) are demonstrated in [23]. In the different direction, 
motivated by the Hermite problem - to represent real cubic irrationals as periodic sequences 
of integers - Beaver and Garrity [1] introduced a two-dimensional analogue of ?(a;). They 
showed that periodicity of Farey iteration corresponds to a class of cubic irrationals, and that 
the two dimensional analogue of ?(a;) possesses similar singularity properties. Nevertheless, 
the problem of Hermite still remains open. Salem in [26] proved (see also [13|) that ?(x) 
satisfies the Lipschitz condition of order where 7 = ^-y^, and this is in fact the best 

possible exponent for Lipschitz condition. The Fourier-Stieltjes coefficients of ?(a;), defined 
as Jq e^™^ d7 (x) , where investigated in [26]. The author, as an application of Wiener's 
theorem about Fourier series, gives average results on these coefficients without giving an 
answer to yet unsolved problem whether these coefficients vanish, as n — > 00. It is worth 
noting that in Section 8 we will encounter analogous coefficients (see Proposition 3). 
Recently, Calkin and Wilf [6] defined a binary tree which is generated by the iteration 

a a a + b 

b a + b ' 6 ' 

starting from the root j. Elementary considerations show that this tree contains every 
positive rational number once and only once, each being represented as a reduced fraction 
[6]. The first four iterations lead to 




Thus, the nth generation consists of 2" ^ positive rationals. It is surprising that the iteration 
discovered by M. Newman [21] 



Xi = 1, Xn+l = l/{2[Xn] + I - Xn), 

produces exactly rationals of this tree, reading them line-by-line, and thus gives an example 
of a simple recurrence which produces all positive rationals. Recently, the authors of [lOj 
produced a natural analogue of this tree, replacing integers r with polynomials r G (Z/2Z)[x]. 
One of the results is that these polynomials also satisfy analogous recurrence (minding 
a proper definition of integral part of rational function, which comes for the Euclidean 
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algorithm). It is important to note that the nth generation of the Calkin- Wilf binary tree 
consists of exactly those rational numbers, whose elements of the continued fraction sum up 
to n. This fact can be easily inherited directly from the definition. First, if rational number 
I is represented as a continued fraction [oq, cti, ctr], then the map | ^ maps | to 
[ao + 1, flr]- Second, the map | ^ maps | to [0, ai + 1, a^] in case f < 1, and to 
[1, oo, ai, Or] in case f > 1- This is an important fact which makes the investigations of 
rational numbers according to their position in the Calkin- Wilf tree highly motivated from 
the perspective of metric number theory and dynamics of continued fractions. The sequence 
of numerators 

0,1,1,2,1,3,2,3,1,4,3,5,2,5,3,4,1,... 

is called the Stern diatomic sequence and was introduced in [28j. It satisfies the recurrence 
relations 

s(0) = 0, s(l) = 1, s{2n) = s{n), s{2n + I) = s{n) + s{n + I) . 

This sequence and the pairs {s{n),s{n + 1)) have also investigated by Reznick [25]. It is 
not surprising (minding the relation to Farey tree) that the "distribution" of numerators, 
which are defined via the moments Q^^^ = Yl^n=2^ +i ^'^'^ i^) i t > 0, has an interesting 
application in thermodynamics and spin physics [H], [7]. 

In [T] it was shown that each generation of the Calkin- Wilf tree possesses a distribution 
function F„(a;), and that -F„(x) converges uniformly to F{x). This is, of course, a well 
known fact about the Farey tree. The function F{x) as a distribution function is uniquely 
determined by the functional equation [1] 

1 F(t^) if 0<x<l. 

This implies F{x) + F{l/x) = 1. The mean value of F{x) has been investigated by several 
authors, and was proved to be 3/2 ([1], [25]). 

On the other hand, almost all the results mentioned reveal the properties of the Minkowski 
question mark function as function itself. Nevertheless, the final goal and motivation of 
papers [1], [2] and this work is to show that in fact there exist several unique and very 
interesting analytic objects associated with F{x) which encode a great deal of essential 
information about it. These objects will be introduced in Section 2. Lastly, and most im- 
portantly, let us point out that, surprisingly, there are striking similarities and analogies 
between the results proved here as well as in [2], with Lewis-Zagier's [19] results on period 
functions for Maass wave forms. That work is an expanded and clarified exposition of an 
earlier paper by Lewis [18]. The concise exposition of these objects, their properties and 
relations to Selberg zeta function can be found in [30j . The reader who is not indifferent to 
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the beauty of Minkowski's question mark function is strongly urged to compare results in 
this work with those in [12]. Thus, instead of making quite numerous references to [19J at 
various stages of the work (mainly in Sections 2, 3, 8 and 9), it is more useful to give a ta- 
ble of most important functions encountered there, juxtaposed with analogous object in this 
work. Here is the summary (the notations on the right will be explained in Sections 2 and 9). 



Maass wave form 


u{z) 




Periodic function on the real line 


Period function 




G{z) 


Dyadic period function 


Distribution 


U{x) dx 


dF{x) 


Minkowski's "question mark" 


L— functions 


Lo{p),Li{p) 


(Mis) 


Dyadic zeta function 


Entire function 


9{w) 


m{t) 


Generating function of moments 


Entire function 


(j){w) 


M{t) 


Generating function of moments 


Spectral parameter 


s 


i- 1 


Analogue of spectral parameter 



As a matter of fact, the first entry is the only one where the analogy is not precise. Indeed, 
the distribution U{x) is the limit value of the Maass wave form u{x + iy) on the real line (as 
y — > +0), in the sense that u{x + iy) ~ y^~^U{x) + y^U{x), whereas "^{x) is the same F{x) 
made periodic. As far as the last entry of the table in concerned, the "analogue" of spec- 
tral parameter, sometimes this role is played by 1, sometimes by |. This occurs, obviously, 
because the relation between Maass form and F{x) is only the analogy which is not strictly 
defined. 

This work is organized as follows. In Section 2 we give a summary of the previous results, 
obtained in p] and [2j. In Section 3 we give a short proof of three term functional equation 
( IT2l) . and prove the existence of certain distributions, which can be thought as close relatives 
of F{x). In Section 4 it is demonstrated that there are linear relations among moments Ml, 
and they are presented in an explicit manner. Moreover, we formulate a conjecture, based 
on the analogy with periods, that these are the only possible relations. In Section 5, the 
estimate for the moments rriL is proved. As a consequence, lim^^^oo = — 2yTog^. In 

Section 6 we prove the exactness of a certain sequence of functional vector spaces and linear 
maps related to F(x) in an essential way. Section 7 is devoted to calculation of number of 
integrals, giving a rare example of Stieltjes integral, involving the question mark function, 
that "can" be calculated. In Section 8 we compute the Fourier expansion of F{x). It is shown 
that this establishes yet another relation among m(t), G{z) and F{x) via Taylor coefficients 
and special values. In the penultimate Section 9, the associated Dirichlet series Ca^(s) is 
introduced. In the last section, some concluding remarks are presented, regarding future 
research; relations between F{x) and the Calkin- Wilf tree (and the Farey tree as well) to the 
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known objects are established. Note also that we use the word "distribution" to describe a 
monotone function on [0, cxd) with variation 1, and also for a continuous linear functional on 
some space of analytic functions. In each case the meaning should be clear from the context. 




Figure 1. Minkowski's question mark 7{x), x e [0, 1] 



2. Summary of previous results 
This section provides a summary of previous results. Let 



OO OO 1 

= y" dF(x) , rriL^ J (^^) dF(x) = 2^2:^ dF(x) . (3) 



1 

Ml 



Both sequences are of definite number-theoretical significance because 



i=lim2i-" V [ao,ai,..,a3]^, = lim 22"" V [0, ai, .., a^]^, (4) 

n-^oo ' * n— »oo ' 

ao+ai+...+as=n ai+...+as=n 

(the summation takes place over rational numbers presented as continued fractions; thus, 
tti > 1 and as > 2). In fact, clarification of their nature was the initial main motivation 
for our work. We define the exponential generating functions M{t) = ^^=o ^^^^ m{t) = 
Er=o TT*^- Thus, 

OO OO 1 

M{t) = J e"*dF(a;), m{t) = J exp (^^) dF{x) = 2 J e^'dF{x). 
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One easily verifies that m{t) is an entire function and tliat M{t) fias radius of convergence 
log 2. There are natural relations among values Ml and rriL, independent of a specific 
distribution, like F{x). They encode the relations among functions x^, L G Nq, and functions 
(;^)^ L e No, given by = E (^i) (tTiY ■ Therefore, 



s>L 



s>L 

On the other hand, the intrinsic information about F{x) is encoded in the relation 

L-l 



„_n \^ / 



(6) 



Further, we have 



M{t) = -^—m{t), m(t) = e^m{-t). (7) 
2 — e 

The first relation is equivalent to the system ([6]), and it encodes all the information about 
F{x) (provided we take into account the natural relations just mentioned). The second one 
represents only the symmetry property, given by F{x)+F{l/x) = 1. One of the main results 
about m(t) is that it is uniquely determined by the regularity condition m{—t) <^ e~^*^°s^, 
as t oo, the boundary condition m(0) = 1, and the integral equation 

oo 

m(-s) = (2e" - 1) y m'(-t) Jo(2v^) dt, s e R+. (8) 



(Here Jq{*) stands for the Bessel function Jq{z) = ^ JJ" cos(z sinx) dx). This equation can 
be rewritten as a second type Fredholm integral equation ([15j, chapter 9). In fact, if we 
denote 

then one has 



oo oo 



m(.) =£(.)- / n^{t)K{s,t)dt, where £(.) = | dt. (9) 







On the other hand, all the results about exponential generating function can be restated in 



oo 



terms of generating function of moments. Let G{z) = ^2 ''^l+iz for \z\ < 1. Moreover, the 

L=0 
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functional equation for Gx{z) implies that there exist all left derivatives of Gx{z) at z = 1. 
Then the integral 

oo 1 



G{z) = I ^^±^ dF{x) = 2 / ^ dF{x). (10) 



x+l 





extends G{z) to the cut plane C \ (1, oo). The generating function of moments does not 
exists due to the factorial growth of M^, but the generating function can still be defined in 
the cut plane C' = C \ (0, oo) by dF{x). In fact, this integral equals to G{z + 1), 

which is the consequence of an algebraic identity 



X 

X 



x+l 



1-xz l--^^{z + l)- 
The following result was proved in [2]. 

Theorem 1. The function G{z), defined initially as a power series, has an analytic contin- 
uation to the cut plane C \ (1, oo) via / fJOj) . It satisfies the functional equation 

^ ^ g(-^)+2G{z + 1) = G{z), (11) 



1-z (l-zy 
and also the symmetry property 

1 „/l \ 1 



G{z + 1) = --G(- + 1 

z^ \z 



z 



Moreover, G{z) ^ 0, if z oo and the distance from z to a half line [0, oo) tends to infinity. 
Conversely, the function having these properties is unique. 

Note that two functional equations for G{z) can be merged into single one. It is easy to 
check that the equation 

- + ^g(-)+2G{z + 1) = G{z) (12) 



z z^ \z 



is equivalent to both of them together. In fact, the change z z'wi the last equation gives 
the symmetry property, and application of it to the term GiXj z^ in ( fT2l) gives the functional 
equation in the Theorem 1. Nevertheless, it is sometimes convenient to separate f|T2|) into 
two equations. The reason for this is that in ffTTj) all arguments belong simultaneously to H, 
M, or H", whereas in (fT2ll they are mixed. This will become crucial later (see the last the 
section) . 

The transition m(t) — > Gi^z^ is given by Laplace transform: 



oo 

\^zG{z) = j m{zt)e~*dt. 
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The same transform applied to the eigen-functions of the Fredholm integral equation ([9]), 
yields the following result [2]. 

Theorem 2. For every eigen-value A of the integral operator associated with the kernel 
K{s,t), there exists at least one holomorphic function Gx (defined for 2; G C \ M>ij, such 
that the following holds: 

2Gxiz + 1) = Gxiz) + ^Gx (i) . (13) 

Moreover, G\{z) for < satisfies all regularity conditions, imposed by it being an image 
under Laplace transform ([17 \. p. 469). 

Conversely: for every X, such that there exists a function, which satisfies i\13^ and these 
conditions, A is the eigen-value of this operator. The set of all possible A 's is countable, and 
A„ — * 0, as n ^ 00. 

Figure 2 shows the functions Gx{z) (for the first six eigen- values) for real z in the interval 
[— 1, —0.2]. The choice of this interval is motivated by Theorem 5. Note also that functional 
equation implies Gx{0) = (| + jx)^>'(~^)- Thus, one has ^^^^j^ — > cxd, as A — > 0. This can 
also be seen empirically from Figure 2. 

Summarizing, there are three objects associated with the Minkowski question mark function. 

• Distribution F{x) = Functional equations ([^ + Continuity. 

• Period function G{z) = Three term functional equation [T^) + Mild growth condition 
(as in Theorem 1). 

• Exponential generating function m{t) = Integral equation ^ + boundary value and 
diminishing condition on the negative real line. 

Each of them is characterized by the functional equation, and subject to some regularity 
conditions, it is unique, and thus arises exactly from F{x). The objects are described via 
the "equality" Function = Equation + Condition. This means that the object on the left 
possesses both features; conversely - any object with these properties is necessarily the 
function on the left. 

As expected, here we encounter the phenomena of "bootstrapping": in all cases, regularity 
conditions can be significantly relaxed, and they are sufficient for the uniqueness, which 
automatically implies stronger regularity conditions. Here we show the rough picture of this 
phenomena. In each case, we suppose that the object satisfies the corresponding functional 
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equation. For the details, see [2]. 

(z) F{x) is continuous at one point =^ F{x) is continuous. 

(ii) There exists e < 1 such that for every z with < 0, we have 

G{z -x) = 0(2""^) as X oo ^ G{z) = 0(1^1"^) as dist{z, M+) oo. 
{in) m\-t) = 0{r^) as t ^ oo ^ < e'"^^^ as t -> oo. 

Corresponding converse results were proved in [2]. For F{x), this was in fact the starting 
point of these investigations, since the distribution of rationals in the Calkin- Wilf tree is a 
certain continuous function satisfying ([2]); thus, it is exactly F{x). The converse result for 
m{t) follows from Fredholm alternative, since all eigen- values of the operator ([9]) are strictly 
less than 1 in absolute value. Finally, the converse theorem for G{z) follows from a technical 
detail in the proof, which is the numerical estimate < — ^^^y^ < 1; as a matter of fact, it 
appears that this is essentially the same argument as in the case of m(t), since this constant 
gives the upper bound for the moduli of eigen-values. 

One of the aims of this paper is to clarify the connections among these three objects, and 
to add the final fourth satellite, associated with F{x). Henceforth, we have the complete list: 

• Zeta function Ca^(s) (see definition (24\ ) below) = Functional equation with symmetry 
s —s i^25\) + Regularity behavior in vertical strips. 

In this case, we do not present a proof of a converse result. Indeed, the converse result 
for G{z) is strongly motivated by its relation to Eisenstein series Gi{z) (see [2] and the 
last Section). In the case of Ca^(s), this question is of small importance, and we rather 
concentrate on the direct result and its consequences. 

3. Three term functional equation, distributions Fx{x) 

In this section, we give another proof of (fT2|) . different from the one presented in [2], 
since it is considerably shorter. For our purposes, it is convenient to work in a slightly 
greater generality. Suppose that A G M has the property that there exists a function -Fa(x), 
X e [0, oo), such that 

m{x+l) = ^dFA(x), dFA(^) = ^dFA(x). (14) 

We omitted the word "continuous" in the description of the function intentionally. For a 
moment, consider F\[x) = F{x) with A = —1. Then F_i(x) is certainly continuous. The 
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reason for introducing A will be apparent later. Thus, let 

oo 


Since F{x) + F{l/x) = 1, we see that for A = — 1 this agrees with the definition fllOl) . This 
converges to an analytic function in the cut plane C \ (1, oo). We have 

1 oo 

2Gx{z + 1) = 2 / dFxix) + 2 [ dF{x) = 

J X — Z J X — z 

1 
oo oo 

2 ! [ —^^ dF,(x + l) = 

J - z Vx+1/ / x + l-z 



oo 

2 



z 



oo 



^ + J-^Gx(-] +Gx{z), where a = / dFxix). 







Xz Xz"^ \z> 

For A = —1 and F„i(x) = -F(a;), this gives Theorem 1. Further, suppose A 7^ —1. Then 



// / Oi 01 

dFx{x) = / dFxix) + / dFxix) = -- — ^a = 0. 

010 
Therefore, the last functional equation reads as 



2Gxiz + l) = ^Gx(^)+Gxiz). 



As a matter of fact, there cannot be any reasonable functions Fxix), satisfying (IHj). Never- 
theless, the last functional equation is identical to (fT3!) . Thus, Theorem 2 gives a description 
of all such possible A. This suggests that we can still find certain distributions Fxix). Fur- 
ther, as it was mentioned, —1 is not an-eigen value of (jH]). Due to the minus sign in front of 
the operator, this is exactly the exceptional eigen-value, which is essential in the Fredholm 
alternative. The above proof (rigid at least in case A = —1), surprisingly, proves that the 
next tautological sentence has a certain point: "—1 is not an eigen-value because it is — 1". 
Indeed, we obtain non-homogeneous part of the three term functional equation only because 
A = — 1, since otherwise a = and the equation is homogenic. 

Distributions Fxix) can indeed be strictly defined, at least in the space of functions, 
which are analytic in the disk D = {z : — ^| < |}, including its boundary. This space 
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is equipped with a topology of uniform convergence, and distribution on this space is any 
continuous hnear functional Indeed, since 

1 

^=1 

define a distribution Fx on the space by {z^, Fx) = L > 1, (1, Fx) = 0, and for any 
analytic function B{z) G C'^, B{z) = Y1'l=o^lz^ , by 

OD 

{B,Fx) = Y.hL{z\Fx). 

L=0 

First, (*, Fx) is certainly a linear functional and is properly defined, since the functional 
equation f|T3l) for Gx{z) imphes that it possesses all left derivatives at z = 1; as a consequence, 
the series YlT=i converges for any p G N (see Theorem 3 for the estimates on 

moments m^^). Second, let Bn{z) = XIlLo ^iT'*'^^' n > 1, converge uniformly to B{z) in the 
circle \z\ < 1. Thus, supi^i^^^ \Bn{z) — B{z)\ = r„ ^ 0. Then by Cauchy formula, 

b^' = — i ^dz. 

^ 2tii J z^+^ 

\z\=l 

This obviously imply that — &l| < ^n, L > 0, and therefore {*,Fx) is continuous, and 
hence it is a distribution. Using the condition dFx{x + 1) = |dFx(x), these distributions 
can be extended to other spaces. Summarizing, we have shown that Minkowski question 
mark function has an infinite sequence of "peers" Fx{x) which are also related to continued 
fraction expansion, in somewhat similar manner. F{x) is the only one among them being 
"non-homogeneous" . 

4. Linear relations among moments Ml 

In this section we clarify the nature of linear relations among moments M^. This was 
mentioned in p^j, but not done in explicit form. Note that gives linear relations among 
moments uil- mL = X]f=o (s) L > 0. These linear relations can be written in 
terms of . Despite the fact that second equality of ([7]) is general phenomena for symmetric 
distribution, in conjunction with ([6]) it gives the essential information about F{x). Therefore, 
let us denote 

fV. 

q(x,t) = (2 - e*)e^* - (2e* - 1)6"^* = y]Q„(a;)-. 

^ — ^ nl 

n=l 

We see that Qn{x) are polynomials with integer coefficients and they are given by 

Q„(x) = 2a;" - (x + 1)" - 2(1 - x)" + (-x)". (15) 
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n 




Qn 


(x) 


n 












Qn 


(x) 


1 




2x - 


-3 


5 




2x^ 


- 15x^ 


+ 10x3 


-30x2 


+ 5x - 


-3 


2 




2x - 


-3 


6 




6x^ 


-45x^ 


+ 20x3 


- 45x2 


+ 6x - 


-3 


3 


2x3 - 


- 9x2 _^ 3^ _ 


-3 


7 


2x^ - 21x^ 


+ 21x^ - 


- 105x^ 


+ 35x3 


-63x2 


+ 7x - 


-3 


4 


4x3 - 


ISx^ + 4x - 


-3 


8 


8x^ - 84x^ 


+ 56x^ - 


- 210x^ 


+ 56x3 


- 84x2 


+ 8x - 


- 3 



Moreover, the following statement holds. 



Proposition 1. Polynomials Q„(x) have the following properties: 

{i) Q2„(a;) G Lq(Qi(x), Q3(x), Q2„,-i(x)), n > 1; 
{ii) degQ2„ = 2n - 1, degQ2n-i = 2n - 1, n > 1; 
iiii) Q2n{x) := ^2«(^) + 3 reciprocal: Q2n{x) = x'^'^^'^Q2n 

X 

oo 

(iv) j Q„(x) dF(x) = 0. 



Naturally, it is property {iv) which makes these polynomials very important in the study 
of Minkowski's ?(x). Here L^^i*) denotes the Q— hnear space spanned by the specified poly- 
nomials. 




Proof, (i) Let qe(x,t) = |(q(x, t) + q(x, — t)), and o^o{x,t) = |(q(x, t) — q(x, —t)). Direct 
calculation shows that, if e* = T, then 

2qe = e^*(3 - T - |) + e'^\3 - ^ - 2T), 2q, = e^*(l - T + |) - e-^*(l -^ + 2T). 
This yields 

°° ^2n T — 1 e* — 1 °° ^2n+l 

5^^^"^^) (2^1)! = = = ^EQ2n+i(x)^^^. 

n=l ^ ' n=0 ^ ' 

The multiplier on the right, = tanh(t/2), is independent of x, and this obviously proves 
part (i). Also, part (ii) follows easily from (ITSl) . 

(Hi) Since Q2n(a;) = ^(3x2" — (x + 1)2" — 2(x — 1)2" + 3), the proof is immediate. 
(iv) In fact, ([7]) gives (2 - e^)M{t) = (2e* - l)M(-t). For real |t| < log 2, we have M{t) = 
e^^ dF(x). This implies 

oo oo 

q(x, t) dF{x) = ^— / Qn{x) dF{x) = 0, for \t\ < log 2, 
"=o 
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and this completes the proof. ■ 

Consequently, there exist linear relations among moments M^. Thus, for example, part {iv) 
(in case n = 1 and n = 3) implies 2Mi — 3 = and 2M3 — 9M2 + 3Mi = 3 respectively. The 
exact values of Ml belong to the class of constants, which can be thought as emerging from 
arithmetic-geometric chaos. This resembles the situation concerning polynomial relations 
among various periods. We will not present the definition of a period, which can be found 
in [TB]. In particular, the authors conjecture (and there is no support for possibility that it 
can be proved wrong) that "if a period has two integral representations, then one can pass 
from one formula to another using only additivity, change of variables, and Newton- Leibniz 
formula, in which all functions and domains of integration are algebraic with coefficients 
in Q". Thus, for example, the conjecture predicts the possibility to prove directly that 
ffx2 ,0 9^1 dxdy = f^-, „ , ^=, without knowing that they both are equal to and 
this indeed can be done. Similarly, returning to the topic of this paper, we believe that any 
finite Q— linear relation among constants Ml can be proved simply by applying the func- 
tional equation of F{x), by means of integration by parts and change of variables. The last 
proposition supports this claim. In other words, we believe that there cannot be any other 
miraculous coincidences regarding the values of Ml- More precisely. 



Conjecture. Suppose, G Q, < /c < L, are rational numbers such that J2k=o "^kMk = 
0. Let£= [^]. Then 

L 

^rkX^ e LQ(Qi(a;),Q3(x),...,Q2£+i(a;)). 

fc=0 

This conjecture, if true, should be difficult to prove. It would imply, for example, that 
Ml for L > 2 are irrational. On the other hand, this conjecture seems to be much more 
natural and approachable, compared to similar conjectures regarding arithmetic nature of 
constants emerging from geometric chaos, e.g. spectral values s for Maass wave forms (say, 
for PSL2(Z)), or those coming from arithmetic chaos, like non-trivial zeros of Riemann's (i^)- 
We cannot give any other evidence, save the last proposition, to support this conjecture. 

5. Estimate for the moments tul 

This chapter deals with the asymptotic estimate for the moments m^. This result was not 
obtained before, and in view of the expression (jlj), it is of certain number-theoretic interest. 
This result should be compared with the asymptotic formula for Ml, see fl2T|) and the 
expression just below it. A priori, as it is implied by the fact that the radius of convergence 
of G{z) at z = is 1, and by ([5]), for every e > and p > 1, one has ^ ^ (1 — e)^, 
as L — s> CX3. More precisely, we have 
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Theorem 3. Let C = 6^^"^'°^^ = 0.18917.... Then the following estimate holds, as L ^ oo: 

Both implied constants are absolute. 

Proof. Fix J e N, and choose an increasing sequence of positive real numbers jij < 1, 
1 < J < J. We wiU soon specify in such a way that iij — ^ uniformly as L ^ oo. An 
estimate for m^, is obtained via the defining integral (recall that F[x) + F{l/x) — 1): 

Ml Mj+l oo 



Indeed, in the first integral, the integrand is bounded by 1. In the middle integrals, we choose 
the largest value of integrand, and change bounds of integration to [0, //j+i]. The same is 
done with the last integral, with bounds changed to [0, oo). Now choose fij — for 
some decreasing sequence of constants Cj. Functional equation for F(x) implies F(x + n) — 
1 — 2~" + 2~^F{x), a; > 0. Thus, 1 — F{x) x 2~^, as a; — oo (the implied constants 
being min and max of the function \&(a;); see Figure 3 and Section 8). Using the identity 
F{x) + F{l/x) — 1, we therefore obtain 

j=l CjVL cjVL 

g-VIciiog2 ^ J2 e-^(^+^^+i^°s2) + e-^^. (16) 
Here we need an elementary lemma. 

Lemma 1. For given J eN, there exists a unique sequence of positive real numbers cl, ...,c*j, 
such that 

* 1 * 1 * 1 * 1 

Cl = -J + C2=- + C3 = ... = -^ + C^ = -. 
<-l <-2 <-J-l <-J 

Moreover, this sequence {c*, 1 < j < J} is decreasing, and it is given by 
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Proof. Indeed, we see that c* = a; determines the sequence c* uniquely. First, C2 = x — ^ = 
Let Fi{x) = X, F2{x) = — 1. Suppose we have shown that cj = p ^ ^^^^-^ for certain 



sequence of polynomials. Then from the above equations one obtains 

xFj{x) — Fj_i{x) 

Thus, using induction we see that Cj = p-^-^^, where polynomials Fj{x) are given by the 
initial values Fq{x) = 1, Fi{x) = x and then for j > 1 recurrently by Fj+i{x) = xFj{x) — 
Fj_i{x). This shows that Fj{2x) = Uj{x), where U{x) stand for the classical Chebyshev 
f/— polynomials, given by 

Uj{cose) = — — — . 

sm b' 

The last equation ct = 4r implies Fj^i(x) = 0. Thus, f/j+i(x/2) = 0, and all possible values 
of cl are given by = a; = 2 cos j^, k = 1, 2, J + 1. Thus, 

. ^ FA^) ^ UA^m ^ sin^ 
' F,M^) f/,-i(x/2) sin|^ ■ 

Since our concern is only positive solutions, this gives the last statement of lemma. Finally, 
monotonicity is easily verifiable. Indeed, system of equations imply < c^, and then we act 
by induction. □ 

Thus, c* > 2 - for some constant 6 > 0. Returning to the proof of the Theorem, for given 
J, let c* be the sequence in lemma, and let c* = CjVlog2. Thus, 

Cl log 2 = — + C2 log 2 = — + C3 log 2 = ... = — h cj log 2 = — . 

Cl C2 Cj_i Cj 

Choosing exactly this sequence for the estimate f|T6l) . and using the bound for cl, we get: 
rriL < (J+ l)e-^=^'°s2 <- {J + l)C^e'^'^ . 

Finally, the choice J = [L^^'^] establishes the upper bound. 
The lower estimate is immediate. In fact, let /x = Then 



1 \^ . / 1 \^ 





The choice c = log~^^^ 2 gives the desired bound. ■ 

The constants in Theorem can also be calculated without great effort, but this is astray of 
the main topic of the paper. 

It should be noted that, if we start directly from the second definition ([3]) of rriL, then 
in the course of the proof of Theorem 3, we use both equalities F{x) + F{l/x) = 1 and 
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2F{^^) = F{x). Since these two determine F{x) uniquely, generally speaking, our estimate 
for itll is characteristic only to F{x). A direct inspection of the proof also reveals that the 
true asymptotic "action" in the second definition ([3]) of takes place in the neighborhood 
of 1. This, obviously, is a general fact for probabilistic distributions with proper support 
on the interval [0,1]. Additionally, calculations show that the sequence ttil/ (L^^'^C^) is 
monotonically increasing. If this is indeed the case, there exists a limit A = lim^^^oo jrjj^^- 
Such a constant A might be of definite interest in metric number theory. This is the main 
topic of the paper [3]. 

As a final remark, we note that the result of Theorem must be considered in conjunction 
with the linear relations rriL = X]s=o (s) (~^)'*'"^'" — O5 the natural inequalities, 
imposed by the fact that is a sequence of moments of probabilistic distribution with 
support on the interval [0, 1]. We thus have Hausdorff conditions, which state that for all 
non-negative integers m and n, one has 

1 



2 f a;"(l - xr dF{x) = ("^^ 



This is, of course, the consequence of monotonicity of F{x). 

6. Exact sequence 

In this section, we prove the exactness of a sequence of continuous linear maps, intricately 
related to Minkowski's question mark function F{x). Let C[0, 1] denote the space of contin- 
uous, complex-valued functions on the interval [0, 1] with supremum norm. For / G C[0, 1], 
one has the identity 

j fix)dF{x) = J2 / /(^)2""dF(x), (17) 


Indeed, using functional equation (l2l), we have 

1 00 00 ^ 

J f{x)dF{x) = I fQdF{x) = J2 I + 
1 

which is exactly f|T7|) . Let denote, as before, the space of analytic functions in the 
disk D = I2; — || < |, including its boundary. We equip this space with the topology 
of uniform convergence (as a matter of fact, we have a wider choice of spaces; this one is 
chosen as an important example). Now, consider a continuous functional on given by 
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T(/) = Jq /(x) dF(x), and a continuous non-compact linear operator [Cf]{x) = f{x) — 
Yl'^=i /(;5;^)2~"'- Finally, let i stand for the natural inclusion i : C C^. 

Theorem 4. The following sequence of maps is exact: 

^ C^C^^ ^ 0. 

Proof. First, i is obviously a monomorphism. Let / G Kei{C). This means that f{x) = 
E^=i/(^)2-". Let xo e [0,1] be such that |/(xo)| = sup \f{x)\. Since ^^=1 2"" = 1, 

this yields /(^^) = /(a^o) for n G N. By induction, /([O, ni, ^2, n/ + xq]) = f{xo) for 
all / G N, and all rij G N, 1 < z < /; here [*] stands for the (regular) continued fraction. 
Since this set is everywhere dense in [0, 1] and / is continuous, this forces f{x) = const for 
X G [0, 1]. Due to the analytic continuation, this is valid for x G D as well. Hence, we have 
the exactness at term (*). 

Next, T is obviously an epimorphism. Further, identity f[T7l) implies that lm{C) C Ker(T). 
The task is to show that indeed we have an equality. At this stage, we need the following 
lemma. Denote [Sf]{x) = J^Zi f d^P"" ■ 

Lemma 2. Let f G C^. Then = 2T(/) + 0(7-2") for x G D. 

Here T{f) stands for the constant function, the golden section, and the hound 

implied by O is uniform for x G D. 

Proof. In fact, lemma is true for any function with continuous derivative. Let x G D. We 
have 

00 

[S'-fKx) = 2-("^+"^+-+'^'-)/([0,ni,n2,...,n, + x]). 

ni,n2,...,nr=l 

The direct inspection of this expression and ([1]) shows that this is exactly twice the Rie- 

1 

mann sum for the integral / /(x) dF(x), corresponding to the division of unit interval into 



intervals with endpoints being [0, ni, n2, n^], G N. From the basic properties of Mobius 
transformations we inherit that the set [0, ni, n2, + x] for x G D is a circle whose 
diagonal is one of these intervals, say J^. For fixed r, the largest of these intervals has end- 
points and where Fr stands for the usual Fibonacci sequence. Thus, its length is 
p p ~ C7~2^. Let Xo,Xi G B^, and sup^j^j^ |/'(x)| = A. We have 

sup |/(xo)-/(xi)| <Ac7-'^ 
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Thus, the Riemann sum deviates from the Riemann integral no more than 

oo 

|[57](x) - 2T(/)| < Ac^-^"" 2-("i+"2+-+-.) = Ac-i-^\ 

ni,n2,...,nr=l 

This proves lemma. □ 

Thus, let / G Ker(T). All we need is to show that the equation f = g — Sg has a solution 
g e C^. Indeed, let ^ = / + Er=i'5"/- % the above lemma, = 0(7"2"). Thus, the 

series defining g converges uniformly and hence g is an analytic function. Finally, g — Sg = /; 
this shows that Ker(T) C lm{C) and exactness at term (**) is proved. ■ 

The eigen-functions of S acting on the space are given by G*{—x) = /g ^ Gx{z) dz + 
jli^xiz) dz (see equations (fT9!) and (l20l) in the next section). Thus, the problem of con- 
vergence of S"'f is completely analogous to the problem of convergence for the iterates 
of Gauss-Kuzmin-Wirsing operator. Let us remind that if / G C[0, 1], it is given by 
[^/](^) = (x+n)'^ f ix+n) ■ Dominant eigen- value 1 correspond to an eigen-function 

j^. As it was proved by Kuzmin, provided that f{x) has a continuous derivative, there 
exists c > 0, such that 

1 

= + O(e-^v^), as n ^ cx); A = [ f{x) dx. 

1 + X log 2 y 



The proof can be found in [;13j. Note that this was already conjectured by Gauss, but he did 
not give the proof nor for the main neither for the error term. For the most important case, 
when f{x) = 1, Levy established the error term of the form 0(C") for C = 0.7. Finally, 
Wirsing [29] gave the exact result in terms of eigen-functions of W, establishing the error 
term of the form {—c)^'${x) + 0{x{l — x)fi^), where c = —0.303663... is the sub-dominant 
eigen- value (Gauss-Kuzmin-Wirsing constant), \E'(x) is a corresponding eigen-function, and 
/i < |c|. Returning to our case, we have completely analogous situation: operator W is 
replaced by S, and the measure dx is replaced by dF{x). The leading eigen- value 1 corre- 
sponds to the constant function. However strange, Wirsing did not notice that eigen-values 
of W are in fact eigen-values of certain Hilbert-Schmidt operator. This was later clarified by 
Bobenko [5] . Recently, Gauss-Kuzmin-Levy theorem was generalized by Manin and MarcoUi 
in [20] ■ The paper is very rich in ideas and results; in particular, in sheds a new light on the 
theorem just mentioned. 

Concerning spaces for which Theorem 4 holds, we can investigate space the C[0, 1] as well. 
However, if / G C[0, 1] and / G Ker(T), the significant difficulty arises in proving uniform 
convergence of the series ^^q'^"/- Moreover, operator S, acting on the space C[0, 1], has 
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ra-1 

additional point spectrum apart from A. Indeed, let Pn{y) = + Yl ^iV^ t)e a polynomial 

2 = 

of degree n, which satisfies yet another variation of three term functional equation 

On 

for certain The comparison of leading terms shows that (5„ = 2^131' ^^^^ indeed for 
this 6n there exists a unique polynomial, since each coefficient aj can be uniquely determined 
with the knowledge of coefficients for i > j. Thus, 

Piiy) = y-h PM=y'-ly + h 

PM = y'- %y' + f,y-i,, pM =y'- §y' + i^^ -Miv + Me- 

The equation for P„ implies that 

SnPniy) = j:^H^)- 
n=l 

Then we have 

Proposition 2. The function Pn{F{x)) is the eig en- function ofS, acting on the space C[0, 1], 
and eigen-value ^i+i^^^ belongs to the point spectrum of S. 

Proof. Indeed, 

|5(P„ o F)l(.) = f lp„ „ f(-I-) i lp„(l - F(. + „)) ^ 

n=l 

lp„ (2- - 2-"F(x)) 5.P„(F(x)).B 

ra=l 

Thus, operator S behaves differently in spaces C[0, 1] and C^. We postpone the analysis of 
this operator in various spaces for the consequent paper. 



7. Integrals, involving F{x) 

In this section we calculate certain integrals. Only rarely it is possible to express an 
integral involving F{x) in closed form. In fact, all results we possess come from the identity 
Ml = |, and any iteration of identities similar to (ITTI) . The following theorem adds identities 
of quite a different sort. 
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Theorem 5. Let Gx{z) be any function, which satisfies the hypotheses of Theorem 1. Then 



1 1 
(^) ^ j Gx{-x)dx = j Gx{-x)F{x)dx- 



11 1 
(,:,:) - / log . dF(.) = 2 / log(l + .) dF(x) = / G(-.) d.; 



1 

2, . 1 



(m) / G{-x){l + x') dF(x) = -; 





1 



/2 
GA(-x)(l-y)dF(x)=0. 



Proof. We first prove identity (i). By f[T^ . for every integer n > 1, we liave 

2Gx{-z-n + l)-Gx{-z-n) = —^——rGx( ^ 

Mz+nY V 



A(2 + ny \ z + n 

Divide tliis by 2^ and, sum over n > 1. By Tlieorem 2, tlie sum on tlie left is absolutely 
convergent. Thus, 



oo 



Gx[-z) = ^ ,^-———,Gx 



Let G\{x) = Gx{z) dz. In terms of G\{x), the last identity reads as 



Gli-x) = y -^Gl ( - - y -^Gl (--). (19) 



n=l n=l 



In particular, setting x = 1, one obtains 

oo ^ 11 



n=l 
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Now we are able to calculate the following integral (we use integration by parts in Stieltjes 
integral twice). 



1 i i 

J Gxi-x)F{x) ^^ = ~ j -^G\{-x)F{x) Ax = -]^G\{-1) + j G\{-x) dF{x 



in 





n=l n=l Q 

1 1 

1 , , . 11 . 1 /■ . . _ , . 1 



^G\-l) + -(- - l)G\{-l) - - / G\{-x) dF{x) = -G^(-l) - - / G,{-x)Fix) dx. 







Thus, the same integral is on the both sides, and 



1 

G,{-x)F{x)dx = -j^GU-1). 







This establishes the statement (z). 

Now we proceed with second identity. Integral (ITU]) and Fubini theorem imply 



1 1 



j G{-z) = 2 y J dF{x) = 2 j log(l + x) dF{x). 



Lastly, we apply (fTTIl twice to obtain the needed equality. Indeed, 

/ = j log(l + X) dF{x) ^Y.Ynj^^^i} + ^) d^(^) = 


— / log(a; + n) dF(x) — I = — loga;dF(x) — /. 
"o 

This finishes the proof of {ii). 

In proving {in), we can be more concise, since the pattern of the proof goes along the same 
line. One has 

oo -| -J oo 

G(-z) = -y ^ ^ G ( - + y — 

n=l ^ ' n=l ^ ' 
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Thus, 



oo 1 



'^=1 

oo ^ 1 1 



Since xdF{x) = ^ = \i this finishes the proof of {iii). Part {iv) is completely analogous. 
■ 

Part {in)., unfortunately, gives no new information about m^. Indeed, the identity can be 
rewritten as 

CO ^ 

which, after regrouping, turns into the identity mQ-mi = |. 

Concerning part {iv), and taking into account Theorem 4, one could expect that in fact 
Ker(T) is equal to the closure of vector space spanned by functions Gx{—x){l — ^). If this 
is the case, then these functions, along with G{z){l + x'^), produce a Schauder basis for C^. 
Thus, if x^ = a^^^Gx{—x){l — y-), then a^"^"* = 2mL. We hope to return to this point in 
a future paper. 

Concerning (i), note that the values of both integrals depend on the normalization of Gx, 

since it is an eigen-function. Replacing GA(2;)by cGx{z) for some c G M, we deduce that 

1 

the left integral is equal to 1 or 0. Then (i) states that J F{x)Gxi—x) dx = or 0. The 



presence of A + 1 in the denominator should come as no surprise, minding that A is the 
eigen-value of Hilbert-Schimdt operator. The Fredholm alternative gives us a way of solving 
the integral equation in terms of eigen-functions. Since |A| < Ai = 0.25553210... < 1, the 
integral equation is a posteriori solvable, and A + 1 appears in the denominators. Curiously, 
it is possible to approach this identity numerically. One of the motivations is to check its 
validity, since the result heavily depends on the validity of almost all the preceding results in 
The left integral causes no problems, since Taylor coefficients of Gx{z) can be obtained 
at high precision as an eigen-vector of a finite matrix, which is the truncation of an infinite 
one. On the other hand, the right integral can be evaluated with less precision, since it 
involves F{x), and thus requires more time and space consuming continued fractions algo- 
rithm. Nevertheless, the author of this paper have checked it with completely satisfactory 
outcome, confirming the validity. 



24 



GIEDRIUS ALKAUSKAS 



Just as interestingly, results {i) and {iv) can be though as a reflection of a "pair-correlation" 
between eigen- values A and eigen- value —1 (see Section 3 for some remarks on this topic). 
Moreover, minding properties of distributions F^{x) (here fi simply means another eigen- 
value), the following result can be obtained. Given the conditions enforced on by (fT^ . 
identity (fTTIl is replaced by (rigid for / G C'^) 
1 ^ 1 

//wdw=-if://(-i-)2-"dF,w. 



Then our trick works smoothly again, and this yields an identity 

1 

j Gx{-x){\ + f^x^)dF^{x) = 0. 



The fact, consequently, is an interesting example of "pair correlation" between eigen- values 
of Hilbert-Schmidt operator in (Q. Using definition of distribution F^, the last identity is 
equivalent to 

oo 
L=l 

and thus is a property of "orthogonality" of G\{z). This expression is symmetric regarding 
fi and A. As could be expected, it is void in case /i = A. As a matter of fact, the proof of the 
above identity is fallacious, since the definition of distributions Fx does not imply properties 
iHM (these simply have no meaning). Nevertheless, numerical calculations show that the 
last identity truly holds. We hope to return to this topic in the future. 

8. Fourier series 

Minkowski question mark function F{x), originally defined for x > by ([1]), can be 
extended naturally to M simply by functional equation F{x + 1) = | + ^F{x). Such an 
extension still is given by the expression ([T]), with the difference that oq can be negative 
integer. Naturally, the second functional equation is not preserved for negative x. Thus, we 
have 

2^'+^(F(x + !)-!) = 2^'(^(x) - 1) for x eR. 

So, 2^{F{x) — 1) is a periodic function, which we will denote by — \E'(x). Figure 3 gives 
the graph of \l/(x) for x G [0,1]. Thus, F{x) = — 2~^\l/(x) + 1. Since F{x) is singular, 
the same is true for \I'(x): it is differentiable almost everywhere, and for these points one 
has \l/'(x) = log 2 ■ \I^(a;). As periodic function, it has an associated Fourier series expan- 

oo 

sion \E'(x) ~ ^ c^e^'^*"^. Since F{x) is real function, c_„ = c^, n G Z. Let for n > 1, 
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Figure 3. \I>(x) 



Cn = On + ibn, and ttQ = y. Here we list initial numerical values for c* = c„(21og2 — Anin) 
(see the next proposition for the reason of this normalization). 

cl = 1.428159, cl = +0.128533 - 0.026840?, = -0.262601 + 0.004128i, 

c1 = -0.521907 + 0.148754i, = -0.140524 - 0.021886i, ^ = +0.198742 - 0.013703i, 
c* = -0.334910 - 0.017869i, c* = +0.285790 + 0.003744^, c* = -0.008479 + 0.024012i 

It is important to note that we do not pose the question about the convergence of this 
Fourier series. For instance, the authors of [26] and [23] give examples of singular monotone 
increasing functions f{x), whose Fourier- Stieltjes coefficients jj^ e^'^*"^ d/(x) do not vanish, 
as n ^ oo. In [26|, the author even investigated f{x) =7{x). In our case, the convergence 
problem is far from clear. Nevertheless, in all cases we substitute — 2~^\I^(a;) instead of 
— 1) under integral. Let, for example, W{x) be a continuous function of at most 

N 

polynomial growth, as a; ^ oo, and let \I'Ar(x) = ^ c„e^'^*"^. Then 

n=-N 

oo 1 

/oo „ 
W^(x)((F(x)-l) + 2-^^;v(a;))dx < J]|Vr(r)|2-"- / |2^(F(a;) - 1) + dx. 

'■=0 

Since 2^{F{x) — 1) G £2[0,1], the last integral tends to 0, as — > oo. As it was said, 
this makes the change of {F{x) — 1) into — 2^^'\&(x) under integral legitimate, and this also 
justifies term-by-term integration. Henceforth, we will omit a step of changing \l/(a;) into 
\I^Ar(x), and taking a limit ^ oo. 
A general formula for the Fourier coefficients is given by 
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Proposition 3. Fourier coefficients Cn are related to special values of exponential generating 
function m{t) through the equality 

m(log2-27rm) _^ 
Cn = — ^ — -7—. — , and Cn = 0[n ). 
2 log 2 — 47rzn 



Proof. We have (note that = |): 

1 1 





1 



1 /^x(iog2-2.in) ^^/^x ^ rni\og2 - 2mn) 

-2mnJ ^ ' 21og2-47rm 



log 2 



The last assertion of proposition is obvious. ■ 

This proposition is good example of intrinsic relations among the three functions F{x), G{z) 
and m{t). Indeed, the moments rriL of F{x) give Taylor coefficients of G{z), which are 
proportional (up to the factorial multipher) to Taylor coefficients of m{t). Finally, special 
values of m{t) on a discrete set of vertical hne produce "Fourier coefficients" of F{x). 
Next proposition describes explicit relations among Fourier coefficients and the moments. 
Additionally, in the course of the proof we obtain the expansion of G{z) for negative real z 
in terms of incomplete gamma integrals. 

Proposition 4. For L > 1, one has 

Ml^LIY- -. (21) 

^(log2-27rzn)^ ^ ^ 

Proof. Let z < be fixed negative real. Then integration by parts gives 

oo oo 

«(^ + 1) = / d(f w - 1) = / (^2-* w d. = 



oo 

n=—oo n=— oo 

where 

oo oo(log2— 27rm) 
V(z)= 1 - p-^{log2-27rm) 1 I / I p"?' d?/ 

/(l--)^' log2-2.in J (l-i3ife^)^' 
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Since by our convention z < 0, the function under integral does not have poles for > 0, 
and Jordan's lemma gives 

oo 

Vn{z) = . . f T. ^i;^e-^dy = 



log 2 — luin 



log 2—2mn ' 



■V[- -), where V^(^)= / ^— e"^ dy. 



log 2 — 2Tcin Vlog2 — 27rm/ J (l—yz) 





The function V{z) is defined for the same values of z as G{z + 1) and therefore is defined in 
the cut plane C \ (0, oo). Consequently, this implies 

G{z +l) = y V ( ) . (22) 

^ ^ ^log2-27rm Vlog2-27rm/ ^ ^ 

The formula is only valid for real 2; < 0. The obtained series converges uniformly, since 
|1 — V; — o^n . I > 1 for n G Z and z < 0. Since 

I log 2—2mn ' — 

00 



this gives us the expansion of G{z + 1) on negative real line in terms of incomplete gamma 
integrals. As noted before, and this can be seen from the function G{z) has all left 
derivatives at 2; = 1. Further, {L — 1)— fold differentiation of V{z) gives 



00 

,L-1 



\/(^-i)(^) = L\J ^_^_^e-^dy l^(^-i)(0) = L!(L - 1)!. 


Comparing with (1221) and (jSj), this gives the desired relation among moments and Fourier 
coefficients, as stated in the proposition. ■ 

It is important to compare this expression with the first equahty of ([7j). Indeed, since m(t) 
is entire, that equality via Cauchy residue formula imphes the result obtained in [1], i.e. 

M,~!^(-!-)'l! (23) 



2 log 2 Vlog2 

It is exactly the leading term in fl2Tl) . corresponding to n = 0. 

9. Associated zeta function 

Recall that for complex c and s, c** is multi- valued complex function, defined as e^^°^'^ = 
gs(iog|c|+iarg{c))_ Heuccforth, we fix the branch of the logarithm by requiring that the value 
of argc for c in the right half plane 9ftc > is in the range (— 7r/2, 7r/2). Thus, if s = a + it, 
and if we denote r„ = log 2 + 27rm, then |r~*| = |r„|~'^e*^''^'"" ~ |r„|~'^e^'^*''^ as n ^ ±00. 
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Minding this convention and the identity fl2Tl) . we introduce the zeta function, associated 
with Minkowski question mark function. 

Definition 1. The dyadic zeta function Cm{s) is defined in the half plane > by the 
series 

C^W = E (|og2-2..n). - 

where Cn are Fourier coefficients of and for each n, (log 2 — 2Tiiny is understood in 

the meaning just described. 

Then we have 

Theorem 6. Ca^(s) has an analytic continuation as an entire function to the whole plane 
C, and satisfies the functional equation 

Cm{s)V{s) = ~Cm{-s)T{-s). (25) 

Further, C,m{L) = for L > 1. CMi-s) has trivial zeros for negative integers: C,m{~L) = 
for L > 1, and = {L — 1)\{—1)^Ml. Additionally, Ca^(s) is real on the real line, 

and thus Cm{s) = Cm{s)- The behavior of Ca^(s) in the vertical strips is given by estimate 

\(M{^ + ^t)\<^t-'^-'/'■e^^'^/' 

uniformly for a < a < b, |t| ^ oo. 

As we will see, these properties are immediate (subject to certain regularity conditions) 
for any distribution f{x) with a symmetry property f{x) + f{l/x) = 1. Nevertheless, 
it is a unique characteristic of F{x) that the corresponding zeta function can be given a 
Dirichlet series expansion, like (Ell). We do not give the proof of the converse result, since 
there is no motivation for this. But empirically, we see that this functional equation is 
equivalent exactly to the symmetry property. Additionally, the presence of a Dirichlet series 
expansion yields a functional equation of the kind /(x + 1) = jf{x) + ^. Generally speaking, 
these two together are unique for F{x). Note also that the functional equation implies that 

oo 

Cx(zt)r(l + it) = J x**dF(a;) is real for real t. Figure 4 shows its graph for 1.5 < t < 90. 



Further calculations support the claim that this function has infinitely many zeros on the 
critical line 3fJs = 0. On the other hand, numerical calculations of the contour integral reveal 
that there exists much more zeros apart from these. 

We need one classical integral. 
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Figure 4. CM(ii)r(l + it) 



Lemma 3. Let A be real number, arctan(^) = e (— |, |); and > 0. Then 



oo 

/ a;^~^e~^ cos(Ax) dx = 



cos(0s)r(s). 



The same is valid with cos replaced by sin on both sides. 
This can be found in any extensive table of gamma integrals or tables of Mellin transforms. 

Proof of Theorem. Let for n > 0, arctan = We will calculate the following integral. 
Let > 0. Then integrating by parts and using lemma, one obtains 



J x'd{F{x) -l) = s J x"-'2-^^(a;) dx = sj^^n J x"-i2-V^^"^ 



dx 





oo 



— / ^2a„ cos(27rna:) — 26„sin(27rnx)^2~^dx = 

oo 

2sr(s) I log 2 + 27rm|~'* cos(0„s) — 6„sin(0„s)j = sT{s) 



(log 2 - 27rin)' 



Note that the function x^ dF{x) is clearly analytic and entire. Thus, sT{s) Cm{s) is an 
entire function, and this proves the first statement of the theorem. Since F{x)+F{l/x) — 1, 
this gives x^ dF{x) — x"^ dF{x), and this, in turn, implies functional equation. All 
other statements follow easily from this, our previous results, and known properties of the 
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CA,(s)r(s + i)| < 



and the last statement of the theorem follows from the Stirling's formula for F— function: 
|r((j + ~ A/27rt'^~^/^e~''l*l/2 uniformly for a < a < 6, as \t\ ^ oo. ■ 

At this stage, we will make some remarks, concerning the analogy and differences with 

oo 

the classical results known for the Riemann zeta function ^(s) = X] Let d{x) denote 

n=l " 

the usual theta function: Q[x) = Yl ^™ > 0- The following table summarizes all 

n.eN 

ingredients, which eventually produce the functional equation both for ({s) and (m{s). 



Function 


as) 


(Mis) 


Dirichlet series expansion 
Functional equation 


Periodicity: 9{x + 2) = 9 


Periodicity: F'{x+ 1) = ^F'{x) 
F'{x) = 



Since F{x) is a singular function, its derivative should be considered as a distribution on 
the real line. For this purpose, it is sufficient to consider a distribution U{x) as a derivative 
of a continuous function V{x), for which the scalar product {U,f), defined for functions 
/ G C°°(M) with compact support, equals to —{V,f') = — J f'{x)V{x) dx. Thus, both 

6{x) and 2^F'{x) are periodic distributions. This guarantees that the appropriate Mellin 
transform can be factored into the product of Dirichlet series and gamma factors. Finally, 
the functional equation for the distribution produces the functional equation for the Mellin 
transform. The difference arises from the fact that for 6{x), the functional equation is sym- 
metry property on the imaginary line, whereas for F'{x) we have the symmetry on the real 
line instead. This explains the unusual fact that in fl2^ we have the summation over the 
discrete set of the vertical line, instead of the summation over integers. 

We will finish by proving another result, which links Ca^(s) to the Mellin transform of 
G{—z + l). This can be done using expansion (122|) . but we rather chose a direct way. Let 

oo 

J G{—z + Vjz'^'^ dz = G*{s). Symmetry property of Theorem 1 implies that G{—z + 1) 



has a simple zero, as z oo along the positive real line. Thus, basic properties of Mellin 
transform imply that G*{s) is defined for < 3fJs < 1. For these values of s, we have the 
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following classical integral: 

1 

\ ^C"" / - xY' da; = r(s)r(l - s) 

io 1 + 2; J 



Sm TTS 



Thus, using ffTOj) . we get 

00 00 _^ 00 00 _^ 00 

G*{s) = j j ^^—dF{x) = j j rr^^'^''^'^'^^^^^ " ihT^ / ^'~'^^(^)- 
00 00 

This holds for < < 1. Due to the analytic continuation, this gives 
Proposition 5. For s G C \ Z, we have an identity 

G\s)= Cm{s-1)T{s)--^. 

sinvrs 

Therefore, G*{s) is a meromorphic function, G*{s + 1) = —G*{—s + 1), and reSs=LG*(s) = 
{—1)^Ml_i. This is, of course, the general property of the Mellin transform, since formally 

M 

G{z + 1) = ^^^0 Mlz^-\ Thus, G{z + 1) ~ E Mlz^-^ in the left neighborhood of z = 0. 

L=0 

10. Concluding remarks 

10.1. Dyadic period functions in H. As noted in [2], one encounters the surprising 
fact that in the upper half plane H, the equation (fTTjl is also satisfied by ^Gi{z), where 

2 

Gi{z) stands for the Eisenstein series, which in HI is given by Fourier expansion ^ — 
Stt^ E^-^ (Ti(ri)e^™^. The reason for this is that Gi{z) transforms under the action of 
PSL2(Z) as (see ^) 

Gi{z + 1) = Gi{z), Gi{-l/z) = z^Gi{z) - 2TTiz. 

Let fo{z) = G{z) — ^Gi{z), where G{z) is the function in Theorem 1. Then for z & M, 
fo{z) satisfies the homogeneous form of the three term functional equation ffTTj) : moreover 
fo{z) is bounded, when '^z — > 00. Thus, if f{z) = fo{z), 



^/(^)+2/(. + !) = /(.). 



Therefore, denote by DPF° the C— linear vector space of solutions of this three term func- 
tional equation, which are holomorphic in EI and are bounded at infinity, and call it the 
space of dyadic period functions in the upper half-plane. Consequently, this space is at least 
one-dimensional. If we abandon the growth condition, then the corresponding space DPF is 
infinite-dimensional. This is already true for periodic solutions. Indeed, if f{z) is a periodic 
solution, then f{z) = -^f{—l/z). Let P{z) G C[2], and suppose that i{z) stands, as usu- 
ally, for the j— invariant. Then any modular function of the form j'{z)P{j{z)) satisfies this 
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equation. Additionally, there are non-periodic solutions, given by /o(-2)P(j(z)). Therefore, 
G{z) surprisingly enters the profound domain of classical modular forms and functions for 
PSL2(Z). Moreover, in the space DPF, one establishes the relation between real quadratic 
irrationals (via G{z), Minkowski question mark function F{x) and continued fraction algo- 
rithm), and imaginary quadratic irrationals (via j— invariant and its special values). Hence, 
it is greatly desirable to give the full description and structure of spaces DPF° and DPF. 

10.2. Where should the true arithmetic zeta function come from? Here we present 
some remarks, concerning the zeta function Cm{s)- This object is natural for the question 
mark function - its Dirichlet coefficients are the Fourier coefficients of F{x), and its special 
values at integers are proportional to the moments Ml. Moreover, its relation to G{z), m{t) 
and F{x) is the same as the role of L— series of Maass wave forms against analogous objects 
[30] . Nevertheless, one expects richer arithmetic object associated with Calkin- Wilf tree, 
since the latter consists or rational numbers, and therefore can be canonically embedded 
into the group of I'deles Aq. The adic distribution of rationals in the n— th generation 
of Calkin- Wilf tree was investigated in p. Surprisingly, Eisenstein series Gi{z) yet again 
manifest, as in case of M (see previous subsection). Nevertheless, there is no direct way of 
normalizing moments of the n—th generation in order for them to converge in the p— adic 
norm. There is an exception. As shown in [1], 

J2 [ao,ai,..,a,] =3-2"-2-^, 

a(,+ai+...+as=n 

and thus we have a convergence only in the 2— adic topology, namely to the value — |. The 
investigation of p— adic values of moments is relevant for the following reason. Let us apply 
F{x) to each rational number in the Calkin- Wilf tree. What we obtain is the following: 




Using (j2]), we deduce that this tree starts from the root |, and then inductively each rational 



2' 

r produces two offsprings: | and | + |- One is therefore led to the following 



Task. Produce a natural algorithm, which takes into account p— adic and real properties 
of the above tree, and generates Riemann zeta function C(s). 
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We emphasize that the choice of ({s) is not accidental. In fact, M— distribution of the 
above tree is a uniform one with support [0,1]. Further, there is a natural algorithm to 
produce "characteristic function of ring of integers of M" (that is, e~'^^^) from the uniform 
distribution via the central limit theorem through the expression 

1 1 

//W.--dx= mnjLld.,.../ dw(^i^^). 

(For clarity, here we take the uniform distribution in the interval [—1, 1]). This formula and 

2 

this explanation and treatment of e~'^^ as "characteristic function of the ring of integer of 
M" is borrowed from [12], p. 7. Further, the operator which is invariant under uniform 
measure has the form [W/](x) = |/(f) + |/(f + |)- Indeed, for every / G C[0, 1], one has 
/q^[Z///](x) ds = Jq /(x)dx. The spectral analysis of U shows that its eigenvalues are 2~"', 
n > 0, with corresponding eigen-functions being Bernoulli polynomials Bn{x) [H]. These, as 
is well known from the time of Euler, are intricately related with ({s). Moreover, the partial 
moments of the above tree can be defined as Yl'i=i i^W') ■ These values are also expressed in 
term of Bernoulli polynomials. As we know, there are famous Kummer congruences among 
Bernoulli numbers, which later led to the introduction of the p— adic zeta function Cp{s). 
Thus, the real distribution of the above tree and its spectral decomposition is deeply related 
to the p-adic properties. This justifies the choice in the task of C{s). 

Therefore, returning to Calkin- Wilf tree, one expects that moments can be p— adically in- 
terpolated, and some natural arithmetic zeta function can be introduced, as a "pre-image" 
of ({s) under map F. 
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